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Abstract 

We study the large scale evolution of a scalar lattice excitation u which 
satisfies a discrete wave-equation in three dimensions, uti'y) = — X^ 7 ' a (l ~ 
7 / )ut(7 / ), where 7,7' 6 are lattice sites. We assume that the dispersion 
relation uj associated to the elastic coupling constants 0(7 — 7') is acoustic, 
i.e., it has a singularity of the type \k\ near the vanishing wave vector, k = 0. 

To derive equations that describe the macroscopic energy transport we 
introduce the Wigner transform and change variables so that the spatial and 
temporal scales are of the order of s. In the continuum limit, which is achieved 
by sending the parameter e to 0, the Wigner transform disintegrates into three 
different limit objects: the transform of the weak limit, the H-measure and the 
Wigner-measure. We demonstrate that these three limit objects satisfy a set 
of decoupled transport equations: a wave-equation for the weak limit of the 
rescaled initial data, a dispersive transport equation for the regular limiting 
Wigner measure, and a geometric optics transport equation for the H-measure 
limit of the initial data concentrating to k = 0. 

A simple consequence of our result is the complete characterization of en- 
ergy transport in harmonic lattices with acoustic dispersion relations. 
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1 Introduction. 



The energy transport by atomistic oscillations in crystalline solids is a central ques- 
tion in solid state physics. To the first order approximation, the oscillations can be 
described by a discrete wave equation 



where u(j) G M is composed out of the displacements of the crystal atoms from their 
equilibrium position, as will be discussed in Sec. 11.11 To analyze physically relevant 
properties of the crystal, such as its thermal conductivity, we first need to understand 
how energy is transported within the crystal via purely harmonic vibrations. Such 
transport properties are determined by the dispersion relation u of the crystal, here 
uj(k) = the "hat" denoting a discrete Fourier transform. If lo is not smooth, 

then depending on the wavelength different types of continuum energy transport 
equations can arise. We follow the basic ideas of Luc Tartar who developed in the 
1980s a mathematical framework that can be used to analyze the weak limits of 
certain nonlinear quantities, the energy density being one of them. 

Our main interest is to characterize the macroscopic evolution of the energy density. 
The starting point of our mathematical analysis is the Wigner transform which can 
be interpreted as a "wavenumber resolved" energy density. Let us leave the details 
for Sec. 12.11 and only summarize the main findings here. The Wigner transform 
W e = W £ [ip] of the field ip G ^(^ 3 ) corresponding to a given normal mode allows 
defining the corresponding energy density, e £ = e e [ip], by the formula 



where e > denotes the "lattice spacing" and x G R 3 is a variable which interpolates 
between the points on the scaled lattice el? '. Here T 3 denotes the 3-torus, and we 
identify T 3 = M 3 /Z 3 . Thanks to the equality f R3 dq J J3 dk W £ (x, k) = H(u t =o,Ut=o) 
we are able to identify e e (x) as the energy density at position x G M d . 

A limit of a sequence (W 2 ^ 2 ]), where e tends to 0, is in general given by a non- 
negative Radon- measure \x G M + . For such limit measures to form a sensible ap- 
proximation for the original dynamical system, the convergence property needs to 
be retained in the time-evolution. In addition, for such a description to be useful, fi t 
should also satisfy an autonomous evolution equation. However, this is typically not 
possible if the initial measure concentrates to the singular set of u, i.e., to the points 
where uo is not smooth. Here we will augment the above Wigner transform scheme to 
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encompass the most common type of singularity encountered in solid state physics: 
the case when uj behaves like \k\ near k = 0. Such modes will occur in general within 
crystal models with short range interactions, and they are particularly important as 
they are responsible for sound propagation in the crystal. With some effort, it is 
likely that our results could be extended to cover any dispersion relation for which 
the singular set consists of isolated points. However, we will not explicitly spell out 
the general result here. 

Our result can be seen as a generalization of the analysis in p] where it is shown 
that fit can be computed from fi by solving a dispersive linear transport equation 
provided that fi has no concentrations at wavenumbers k where the dispersion relation 
u is not C 1 . 

If the dispersion relation is not almost everywhere C 1 (with respect to the initial 
Wigner-measure) , we have to resolve finer details of the asymptotic behavior of the 
sequence of initial conditions. More precisely, we show that if the sequence of initial 
excitations is bounded and tight in ^(f^ 3 ) then for all t G M there are two measures, 
a Wigner-measure [it on M. 3 x T 3 and an H- measure /ij 1 on R 3 x S 2 , and an L 2 - 
function <p t such that W £ [ip e (t/e)} converges along a subsequence to (fi t , fif, 4>t) in a 
certain weak sense (Theorem 13.21) . The subsequence can be chosen independently of 
t, and it will only be relevant for determining the limit of the initial data, that is, 
(/io, /ig , 4>o). For all other times i 6 R, the measures ji t , [if and the L 2 -function <p t 
can be determined using the transport equations 

d t fit{x,k) + ^Vuj{k)-V x ^{x,k) = 0, A; Git, (1.2) 
dt$(x, q) + ^V^o(g) • V x tf(x, q) = 0, q G , (1.3) 
d 2 = div(^A o V0), (1.4) 

together with the initial conditions 

v\ t =o = Vo, V U \ t=0 = Vo, 0l i=o = < A), d t 4>(t,q)\ t =o = -iu (q)4>o(q). (1.5) 

Here the constant matrix A and the function u Q are determined by the Hessian of 
the square of the dispersion relation u (defined in (12. 6p ) at k = 0, explicitly 

A = ^DV(O), u (q) = ^q ■ A q. (1.6) 

Equation (11 . 2[) describes the propagation of energy along the harmonic lattice with 
the group velocity Vo>(/c)/(27r). Equation (jl.4p is the wave equation which describes 
the evolution of macroscopic fluctuations. Equation (11.31) is usually known under the 
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name "geometric optics" and it describes the evolution of macroscopic fluctuations 
whose wavelength is much longer than the lattice spacing e and much smaller than 
1, the wavelength of the fluctuations resolved by 4>. 

It follows easily from our analysis that the sum of the energies of /x, /x H and is 
a constant of motion and equals the limiting value of the total energy of the initial 
excitations. 

Standard results in semi-classical analysis show that the Wigner measure fi and 
f R3 dq W 1 [0] (x, q) are non-negative. Hence, a key step in our analysis is to demon- 
strate that also /i H has a sign. The measure /x H is closely related to the H-measure 
which was introduced by L. Tartar in [2] and P. Gerard in [3], but in general it differs 
from the H-measure. 

The Wigner transform, or the Wigner function, was originally introduced to study 
semi-classical behavior in quantum mechanics but it has been proven to be a useful 
tool in studying large scale behavior of wave equations, as well [H [5] . In particular, 
the method of calculating continuous, macroscopic energy by finding the limit object 
of a sequence of energies e e on rescaled lattice models is one that has been widely 
used and justified in, for example, [TJ [6j [7J. In [I], the Wigner transform of the 
normal modes is employed in solving the macroscopic transport of energy in the 
above harmonic systems for deterministic initial data. The same system is considered 
in [8] with random initial data and in a larger function space, however, excluding the 
type of concentration effects we study here. 

The main use of the Wigner transform is that, unlike the energy density e e itself, it 
contains enough information so that it satisfies a closed equation in the limit e — > 0. 
Indeed, it was shown in pQ that, as long as there is no concentration on the singular 
set of the dispersion relation faster than e 1 ^ 2 , the Wigner transform of the time- 
evolved state vector ip(t/e) converges to a limit measure jit on R x K. Here K is a 
suitable compactification of T 3 \ §, § = the singular set, which allows a continuous 
extension of the group velocity Vuj. The measure jit is then proven to satisfy the 
transport equation 

dtjkix, k) + ^Vu(k) ■ V x fl t (x, k) = 0. (1.7) 

It was also shown in [1] , that if the weak limit of the rescaled initial data exists, then 
the limit satisfies the continuum wave equation. However, removing the assumption 
about the rate of concentration, as well as combining the result with a non-zero weak 
limit for the initial data remain open questions. 

In this paper we will show how to overcome the above difficulties in a physically 
relevant class of models with a singular dispersion relation. In Section [2] we will first 
present the microscopic dynamical model in detail. In Section 12.11 we will define the 
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Wigner transform, and discuss its relation to the energy of the microscopic lattice 
model. The main results will be presented in Section [3l 

1.1 Relation with solid state physics 

A crystal in solid state physics is a state of matter in which the atoms retain a nearly 
perfect periodic structure over macroscopic times. The Hamiltonian model used for 
the time-evolution in such a crystal is, to the first order accuracy, harmonic. If we 
assume that each periodic cell of the idealized perfectly periodic crystal structure 
contains n atoms, then we can form a vector 5(7) e M 3n out of the displacements 
of the atoms in the periodic cell labeled by 7 G Z 3 . The (classical) Hamiltonian 
equations of motion of this harmonic model are then 

ft (7. *) = — Pi(7» Pi(l, t) = - V] A(7 - i)i,vqi{l\ t) (1.8) 

mi tji 

where 7 6Z 3 ,i = l,..., 3n, and rrii denotes the mass of the atom whose displacement 
qi measures. 

By the change of variables to ft (7) = mfqi^), ^(7) = m i 2 ^(7), these equations 
can be transformed into a standard form whose force matrix is given by Ai^y)^ = 
m ! ~ 1 ^ 2 A(7)j i! /m ! 7 1 ^. The standard form equations can then be solved by Fourier 
transform, and a diagonalization of the remaining multiplicative evolution equations 
decomposes the 3n vector degrees of freedom into independent normal modes, called 
phonons in solid state physics. Each normal mode is a complex scalar field on 
the crystal lattice, and its time-evolution is unitary and uniquely determined by 
the corresponding dispersion relation Ui(k) on T 3 . More details about the related 
mathematical issues can be found in [TJ [8] . 

From the physical perspective it is highly relevant to understand energy transport 
in the case where u is acoustic in the sense that it contains a \k\ singularity. In 
this case u is degenerate for large wave-numbers in the sense that the velocity at 
which long waves propagate converges to the speed of sound. In other words, the en- 
ergy density can travel ballistically over large distance without experiencing dilution 
effects due to dispersion. 

Acoustic dispersion relations arise in general from atomistic Hamiltonians with 
short range harmonic interactions. For instance, for the type of interactions consid- 
ered in [1] all dispersion relations are Lipschitz continuous and piecewise analytic. 
In solid state physics, the modes are accordingly divided into optical and acoustic 
depending on this regularity: if the dispersion relation is regular at k = 0, then 
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the mode is called optical, if it behaves as \k\ at k = 0, then the mode is called 
acoustic. The latter name arises as these modes are believed to be responsible for 
the propagation of sound waves in the crystal. 

The simplicity of the Fourier-picture is the reason that the dispersion relation is 
the starting point of most of the more physically oriented publications. Obviously, 
within the context of harmonic defect-free crystals both approaches are completely 
equivalent. In the beginning of Section [2] we will provide the mathematical details 
which establish this equivalence on an elementary level. 

Finally, let us remark that the discrete linear wave equation alone does not suffice 
to determine the physically relevant properties of the crystal, such as its thermal 
conductivity. However, it forms the basis for perturbative treatments using which 
these questions can be addressed. We refer to [9j [TOj CD] for further details on the 
physical aspects of the topic and to [TJ] for a review about related open problems. 

Acknowledgments 

We would like to thank Alexander Mielke and Herbert Spohn for several instruc- 
tive discussions. JL was supported by the Deutsche Forschungsgemeinschaft (DFG) 
projects SP 181/19-1 and SP 181/19-2. 

2 The microscopic model 

There are two mathematically equivalent descriptions of harmonic crystals. On the 
one hand, one can work with the Hamiltonian equations of motion and analyze the 
properties of the solutions. At least formally anharmonic crystals can be discussed 
in the same way. On the other hand, one can utilize the periodicity and linearity to 
condense the Hamiltonian into the dispersion relation. While this approach is very 
elegant, it cannot be used to directly analyze nonlinear models. 

We will show in this chapter how for harmonic lattices the first approach reduces 
to the second one. Then we demonstrate how the tool of Wigner functions can be 
brought to bear on the reduced description. 

We assume that the scalar excitation 1^(7), 7 G Z 3 satisfies the discrete wave 
equation 



with initial data (u\ t=0 , d t u\ t=0 ) G X = £ 2 x (-2- The numbers a (7 — 7') are the 




(2.1) 
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Figure 1: Values of 74(7, t) along the axis 72 = 73 = for t = 0, t = 0.1/5, t = 0.9/e 
with e = j * 1CT 1 . The evolution is given by (12.11) with the nearest neighbor elastic 
couplings (I2.12p , and the initial conditions are u t =o(0) = 1, 1^=0(7) = for all 7 7^ 
and u t =o = 0. 



elastic coupling constants between the sites 7 and 7'. We assume that a is real and 
symmetric (a (—7) = a (7)). Clearly system (12.11) can be written in a Hamiltonian 
form and the energy 



E(u) = H(u,ii) 



\ 7 ez 3 7GZ 3 



u(7)ot(7 - 7 / )tt(7 / ) 



(2.2) 



is constant along solutions. Depending on the initial conditions the solutions of 
system ( 12.11) may develop large scale oscillations which carry a finite amount of 
energy, cf. Fig. [1] where snapshots of u at several times are plotted. 

Since system (12.10 is linear and invariant under discrete translations, we can write 
the solutions in a closed form using the Fourier transform. 



Definition 2.1 We define the Fourier transform ^(Z 3 ) - ► L 2 (T 3 ) by extending 



(^*V0(*) = i>( k ) = E e- a "* 7 ^(7) 



(2.3) 



/rom t/> with finite support to all of^i^ 3 )- -Here T 3 = IR 3 /Z 3 denotes the unit 3-torus. 
The inverse transform is pointwise convergently defined by the integral 



{Fk~$){i)= / dke 2 ^(k)=^). 



(2.4) 



T 3 



where the measure dk is induced by the Lebesgue measure on [— ~, ^] 3 . In particular, 



J T 3 dk = 1. 
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If one applies the Fourier transform to equation ( 12.11) one obtains the simpler system 

d_fu(k,t)\_f l\fu(k,t)\ 
dt\v(k,t)J \-co 2 (k) 0j\v(k,t)J- 

The function u : T 3 — > R is the dispersion relation and is related to the atomistic 
Hamiltonian via the following formula 



u(k) = y/a(k) = a(7)cos(2vT7 ■ k). (2.6) 

Y 7 6Z 3 

where we have employed the assumption a(j) = a (—7). Since a is real and satisfies 
the above symmetry property, we find that u is also real and symmetric, i.e. u(k) = 
uo(-k). 

Now diagonalizing the matrix on the right hand side of (12. 5p motivates combining 
the real scalar fields u, v into the two complex fields ip± = ip±[u, v] G ^i^ 3 , C) defined 
by the formula 

$a(k) = -4= (u(k)u(k) + iav(k)) (2.7) 
v2 

where a = ±1. For all (u, v) G X, we clearly have G L 2 (T 3 ), and thus ift a G ^(Z 3 )- 
In addition, since we assumed u>(— k) = u(k), we also have ^-(t) = i J +(l) f° r au 7- 
The transformation can always be inverted by applying 

v = u=-^=(^ + + ^). (2.8) 

These fields are normal modes of the harmonic system, since (12.51) implies that 
"0±(7,£) = 1 P±[ u if)i v satisfy the evolution equations 

d_ ($+(k,t)\ = (u(k) \ (Mk,t)\ 

which are readily solved to yield for all t G R, A; G R 

Vi ± (fc,t) =e Tiw(fc) ^ ± (fc,0). (2.10) 

These are exactly the two evolution equations corresponding to a "phonon" mode 
with a dispersion relation to. 

After these reduction steps it is obvious that the dispersion relation uj fully deter- 
mines the properties of the solutions. We will assume throughout this paper that uj 
is of acoustic type in the following precise sense: 
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Definition 2.2 We call u G C(T 3 , [0, oo)) an acoustic dispersion relation if X = oj 2 

satisfies: 



1. AgC( 3 »(T 3 ,[0,oo)). 

2. A(0) = 0, and the Hessian of X is invertible at 0. 

A dispersion relation is called regular acoustic, if it is acoustic and X(k) > for 
k ^ 0. The 3 x 3-matrix Aq is the Hessian of \X at k = and ujq{o) = y/q ■ A^q. 

These assumptions are fairly general: as discussed in the introduction, all stable 
harmonic interactions have non-negative eigenvalue functions A, and for interactions 
of the type discussed in p] u = VX is Lipschitz continuous. 

A prototype for the kind of dispersion relations we will consider here is the disper- 
sion relation of the nearest neighbor square lattice: 



UQQ (k) = [^2(1 - cos^tt/O) 2 . (2.11 



This is clearly a regular acoustic dispersion relation, in the sense of Definition 
and for it A is proportional to a unit matrix and uj (q) = 27t\q\. The corresponding 
elastic couplings are given by a nn {^' — 7) = — Ay 7 , where A is the discrete Laplacian 
of the square lattice. Explicitly, 

f6, if 7 = 0, 

aw(7) = < -1, if|7l = l, (2-12) 
[0, otherwise. 

To allow the creation of macroscopic oscillations we work with sequences of initial 
conditions that depend on the scaling-parameter e > and consider the asymptotic 
behavior of the solutions as e tends to 0. 



2.1 Energy density and the lattice Wigner transform 

From now we will focus on analyzing asymptotic behavior of the fields ip± as e tends 
to 0. As is carefully discussed in p], generalizing the definitions of the energy density 
and of the Wigner transform to the discrete setting is not completely obvious. In an 
attempt to minimize unnecessary repetition of certain basic results related to Wigner 
transforms, we will resort here to the definitions used in [7] which will allow us to 
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rely on the properties proven in Appendix B of that reference. However, we wish to 
keep in mind that this choice might not be optimal for all purposes, and we refer 
the interested reader to the discussion and to the references in [U [13] for further 
possibilities. 

We employ here the definition that for any state (u,v), its energy density, e e = 
e £ [u, v], scaled to a lattice spacing e > 0, is the tempered distribution defined via the 
complex fields ip a = ip a [u,v) in ( 12. 71) : 

e £ (x)= £ h/v(7)| 2 (2.13) 

7GZ 3 <7=±1 

where 5 denotes the Dirac delta-distribution. This is a manifestly positive distribu- 
tion, and identifiable with a measure whose total mass equals the total energy: 



j dxe £ (x) = \ E \Mi)\ 3 = \ E ii^n 2 

J 7 gz 3 cr=±l <T=±1 

= ^E [ d k\u{k)u{k) +iav(k)\ 2 = H(u,u) < 00. (2.14) 

<T=±1 ^ 

This justifies calling e e an energy density: it defines a distribution of the positive total 
energy between the lattice sites. The symmetry of to implies that |-0-(7)| = I^+Ct)! 
and thus we can also identify the energy density directly with the norm-density of 
^ + : 

e £ (x) = ^5(x-e 7 )l^ + (7)| 2 - (2.15) 

Then also for all e > 0, 

H(u,v) = J dxe £ (x) = 11-0+111(23) = ll^+lli^frs) = ll^-lliafrsj (2.16) 

which is conserved when (u,v) = (u(t),v(t)). 

We have given an example of the time-evolution of the so defined energy density 
in Fig. [21 The last panel in the figure contains the most obvious features which 
are implied by the corresponding macroscopic evolution equation: the points of dis- 
continuity of the solution. The macroscopic initial data is given by fi (dx,dk) = 
S(x)^\u(k)\ 2 dxdk, which has no concentration at k — 0. Thus only ( 11.21) is relevant. 
It is readily solved to yield as the energy density 



e(x, t) 



f dkS(x-t^Vcu(k)) = t- 3 [ dk6{Z-±Vu{k)) (2.17) 
Jt 3 7t 3 
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Figure 2: First three panels: Snapshots of the energy density t) | 2 in the plane 



73 = at t = 0.1/e, t = 0.3/e and t = 0.95/e, e 



128' 



with initial data and elastic 



constants as in Fig. [TJ The plot is of the logarithm of the density, with all values 
less than a fixed cut-off shown white. Last panel: Plot of the restriction to the plane 
x 3 = of the singular set of the solution to the transport equation (jl.2p with the 
corresponding macroscopic initial data. As the solution is scale-invariant, no explicit 
length scale has been denoted. 
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Evaluating such integrals has been considered, for instance, in Sec. 6.4 of pQ. | Va>(fc)| 
has its maximum near the point of discontinuity of the gradient, at k — 0. This 
defines the outer circle outside which the solution must be zero. Inside the circle, 
the solution has a finite density, apart from points which correspond to values of k 
for which the Hessian of to is not invertible. We have computed the positions of such 
points using Mathematica, and plotted the result in the last panel in Fig. [2 For a 
reader interested in the details of the computation, we point out that considering the 
case Xs = simplifies the problem, as it implies that either k$ = or k% = |. 

We are interested in the limiting behavior of e £ [u(t/e),u(t/e)}, as e tends to 0. 
Since the velocity of the waves with wave vector k depends on k it is necessary to 
work with an object that encodes the density of waves with wave vector k G T 3 at 
x G M 3 . This job is conveniently done by the Wigner-transform. In order to avoid 
certain technical difficulties we are going to define our Wigner-transform only in the 
sense of distributions, i.e., via a duality principle. 

First we introduce the space of Schwartz functions. 

Definition 2.3 Let Sd = S(M. d ) denote the Schwartz space, and || • \\s d ,N the corre- 
sponding N :th Schwartz norm. Explicitly, with a denoting an arbitrary multi-index 
and with (x) = yl + x 2 , then 

\\f\\ Sd , N =supmax\(x) N d a f(x)\. (2.18) 
We also employ the shorthand notation S = S%. 

To extract the relevant weak limits from the sequence ip E as e tends to we have to 
specify a space of suitable test functions. Since we want to track the evolution of 
three different kinds of lattice vibrations (short-, medium- and long- wavelength) we 
require a somewhat involved and non-standard notion of multiscale test functions. 

Definition 2.4 We call a test function a G C^(M 3 x T 3 x R 3 ) admissible, if it 

satisfies the following properties: 

!■ su Pk, q ,\a\<N \\d a a(-,k,q)\\ StN < oo, for all N > 0, 

2. q i— > a(x, k, q) is constant for all \k\oo > \ and x G M 3 . 

3. There is a function b G C(°°)(M 3 x T 3 x S 2 ) such that for any N > 

sup \\a(-,k,q) — b(-,k,A)\\ s>N —> 0, when R — > oo. (2.19) 
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The first condition can be summarized as follows: we assume the test-functions to be 
Schwartz in x and smooth with bounded derivatives in k and q. The above require- 
ments are not minimal. The second condition is only needed in order to guarantee 
that k i — > a(x, k, k/e) would always be smooth on T 3 . Also, taking arbitrarily large 
N in the last step is not necessary, most likely iV = d + 3 = 6 would suffice. 

Having the notion of admissible test functions at our disposal we can define the 
central object of this paper: the Wigner transform. 

Definition 2.5 Let ip G ^(Zi 3 )- We define the lattice Wigner transform W £ [ip] at 
scale e > by 

(a,W £ [i;])= [ dp [ dfca(p,fc,*)^(fc-e§Mfc + e§), (2.20) 

JR 3 Jt 3 

where a is an admissible test function and a = T x ^ v a, i.e., 

a(p,k,q)= I dxe- 2lTipx a(x,k,q). (2.21) 
Jr 3 

The L 2 -Wigner transform W^l^cp] of a function <p £ L 2 (M 3 ) at the scale e > is 
given by the distribution 

b ~ (b, W^M) = [ dpdqb(p,q)4>(q-e*)j>(q + ef ) (2.22) 

JR 3 xR 3 

for all b G S(R 3 , C^(R 3 )), and with b = T x ^ v b. 

The test-function space iS(M 3 , C^°°^(1R 3 )) used above to define W COXiX , is obtained via 
the family of seminorms Pat(6) = supi a i i^i i g | <iV \(x) N D a b(x, q)\ with b G C^°°^(1R 3 x 
M 3 ). This is a Frechet-space, and W con t[0] is a continuous functional on it for any 
G L 2 (M 3 ), as the following estimate reveals: 

\(b,W<$M)\ < sn P \(pfb(p,q)\ ||0f / dp(p)-\ (2.23) 
p,q Jr 3 

Although S e is not dense in this test-function space, it is nevertheless enough to know 
how W cont acts on it. More precisely, if Wi = W^lfa], i — 1,2, and (b, W\) = (b, W 2 ) 
for all b G Sq, then W\ = W 2 - This follows straightforwardly from an estimate similar 
to (12.231) using smooth cutoff functions to cut out the infinity of the g-variable. In 
addition, we will also need the property that, if b(x, q) = f(x), with / G S3, then 

(b^W^M = [ dxW)\m\ 2 - (2-24) 
Jr 3 
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That is, at least formally, L 3 dqW^ nt [(f)](x, q) = \<j)(x)\ 2 . 

In [7J the Wigner transform of a lattice state was defined as a distribution Wf att G 
iS'(1R 3 x T 3 ). The above definition is simply a refinement of this definition: formally 

for any ijj 

W £ (x, k,q) = 5(q- \ )Wl tt (x, k). (2.25) 

This follows immediately from Eq. (B.6) of [7J, after one realizes that if a is an 
admissible test function, then (x, k) i— > a(x,k,k/e) belongs to S(M. 3 x T 3 ) for any 
e > 0. This identification immediately allows us to use the results in [7J and to prove 
that many of the basic properties of the usual Wigner transform carry over to the 
multi-scale Wigner transform. Particularly important for us is the following relation: 

Proposition 2.6 For any f G iS(M 3 ) ; the test function af(x,k,q) = f(x) is admis- 
sible, and for all e > and for all {u, v) G X , 

(f,e £ ) = (a f ,W £ [i;}) (2.26) 

where e £ = e £ [u, v] and ip = ip+[u, v]. 

Proof: By inspection, we find that a is a well-defined test-function, and ip G £2- Then, 
by the above mentioned relation, {aj, W e [ip]) = (Jf, Wf att [^]) with Jf(x,k) = f(x). 
On the other hand, it follows directly from the definition of W{ att (equation (B.2) in 
[7J) that 

(^/,^t[#= E V>(7)#7)/ dfce 2 ^^-^/(5(Y + 7)/2) 

7,7 

= E^)i^)i 2 = ^ e£ )- ( 2 - 27 ) 

This proves fl2T26l) . □ 
The above Proposition can be formally summarized by the formula 

e e (x) = [ dk [ dqW £ [ip](x,k,q) (2.28) 

which implies also (in the sense of choosing any suitable test-function sequence ap- 
proaching pointwise 1) 

I dxdkdqW £ {iP](x,k,q) = ||V||| 2 ( 2 3)- (2.29) 

</R 3 xT 3 xR 3 
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As noted earlier, analogous results hold for the L 2 - Wigner transform. Let us also 
remark here that, if the sequence W^ tt [y> e ](-, ■) converges, then the limit is given by a 
non-negative Radon- measure [Tj. The results of the next section will show that also 
P^iitJV^K'; ■ /e) has such a limit property. 



The macroscopic evolution is obtained by sending e to 0. Our objective is to char- 
acterize the asymptotic behavior of the Wigner function W /£ [ , £ ]. The limit strongly 
depends on the dispersion relation u. We will consider here regular acoustic disper- 
sion relations, keeping in mind that in the case of a scalar field and nearest neighbor 
interactions in Z 3 the dispersion relation u is given by (12. lip which is regular acoustic 
with cuo(q) = 27r|g|. The main achievement of this paper is that complicated assump- 
tions concerning the concentrations of the Wigner transform W e in wave-number 
space as e tends to are no longer needed. The only remaining requirements are 
boundedness and tightness of the sequence of initial excitations. 

Assumption 3.1 We consider a sequence of values e > such that e — > 0. For each 
e in the sequence we assume that there is given an initial data vector ipQ G ^(Z 3 ) 
such that 

1. sup £ ||^q|| < oo. 

2. The sequence ipQ is tight on the scale e~ x : 



After these preparations we are in a position to state our result. The main point 
is that if u is a regular acoustic dispersion relation the asymptotic behavior of the 
energy density is characterized by precisely three different objects: the weak limit 
(macroscopic waves), the H-measure (short macroscopic waves) and the Wigner mea- 
sure (microscopic waves) and no assumption concerning energy concentrations except 
those stated in Assumption 13.11 are required. 

Theorem 3.2 Let i}jq G ^(Z 3 ) be a sequence which satisfies Assumption ^ . 1[ Let u 
be a regular acoustic dispersion relation and define tpf G ^(^ 3 ) f or all t E M. by the 
formula 



3 Main results 



lim limsup 





7\>R/e 



r t (k) = e 



(3.2) 
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Let also T 3 = T 3 \ {0}. Then there are positive, bounded Radon measures fJ>o, f$ on 
M 3 x T 3 and R 3 x S 2 , respectively, a function <fto £ L 2 (IR 3 ) and a subsequence (not 
relabeled) such that for all admissible test functions a and i el, 



lim(a, W^ t/£ }) = / dfnix, k) b(x, k, ||t) (3.3) 

£ ^ U ./B3yT3 1 



'XT 



+ / dfxf(x, q) b{x, 0, g) + (a 0) W^tW) > 

where b(x, k, q) = Hindoo a(x, k, Rq) for \q\ = 1, ao(x, q) = a(x, 0, q), and (fit, IH an d 
/if are given by 



(?) := e-^Mq), (3-4) 



/ a(x, k) dfj, t (x, k) := / a(x + £^Vu;(fc), k) dpLo(x, k), (3.5) 
/ 6(x, g) dfif {x, q) := / b{x + t±Vu {q), q) d$(x, q) . (3.6) 



Moreover, for all t the energy equality 



limllVol 2 = M^ 3 x T 3 ) + ^ H (M 3 x S 2 ) + INIiw (3.7) 



holds. 



Remark 3.3 R is immediate from the definition that <fi, /x and /i are u>eafc solutions 
of the set of decoupled linear transport equations ( fLifrfMp . 

The proof of Theorem 13 . 21 also shows that the subsequences which are extracted in the 
statement of the theorem can be characterized by a simple condition. In particular, 
the initial state of the wave-equation, (fio, is determined as the weak-L 2 (IR 3 ) limit of 
the sequence of the functions (<$j) with Fourier-transforms 

m = Sf*i^ ifWU<i, (38) 

I (J otherwise. 

The exact characterization is contained in the following Corollary whose proof will 
be given in Section [5j 
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Corollary 3.4 Let (/0o) be a sequence which satisfies Assumption \3.1[ Suppose that 
0q converges weakly to (fio, and that lim e ^o( a ; W 6 ^]) exists for every admissible 
testfunction a. Then there are unique positive, bounded Radon measures Ho, L^o on 
M 3 x T 3 and IR 3 x S 2 , respectively, such that for every admissible testfunction a 



lim(a, W e [r }} = I dfio(x, k) b{x, k, ±) (3.9) 



+ / d${x, q) b{x, 0, q) + (ao, W^M) , 

where a and b are defined as in Theorem \3.Sl In addition, then 113. 3\) holds for all 
t 6 K along the original sequence e with the initial macroscopic data determined by 
the triplet (/z , $ , 4>o)- 

The measure /i H is closely related to H-measures which have been introduced in the 
context of oscillatory solutions of partial differential equations by L. Tartar [2] and 
P. Gerard [3]. The precise nature of the connection of between /z H and H-measures 
is irrelevant for the purpose of this paper, but for the convenience of the reader we 
include a brief discussion. 

Definition 3.5 (H-measures) Let (f) £ G L 2 (R d ) be a sequence which converges 
weakly to as e — > and let v G .M + (IR d x S^ 1 ) be a nonnegative Radon mea- 
sure. If 



lim / J' :r ^g(ai0 e )J r x .^ (? (a20 £ )V'( A) dk = 

£ ~* JR d 7m° 



dv(x, q) a 1 (x)a 2 (x)ip(q) 
S d-i 



for all ai,a2 G C£°(R ) and ip G C(S ), then v is the H-measure generated by the 
sequence <p £ . 

The connection between /x H and H-measures is established by the following 

Proposition 3.6 Let <p £ G L 2 (IR d ) be a tight sequence converging weakly to as 
e — > oo such that 

lim lim / dA; |0 £ (A;)| 2 = 0. (3.10) 

p _0 e-»0 J £ \ k \> p 

If <p £ generates an H-measure v as e —> and a is an admissible testfunction, then 
lim dx dqa(p, eq, q) (f> £ (q — f)0 E (<? + f) = / dis(x, q) a(x, 0, q). (3.11) 
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Proof: See [H] 



□ 



Note that we will never assume that <p e converges weakly to 0, hence we cannot define 
the H-measure associated to <p e . Moreover, unlike the H-measure v the measure ^ H 
does not take the contribution of oscillations with wavelength 1/e into account. For 
these reasons /x H differs from v in general. 



4 Proof of Theorem 3.2 



Let a be an admissible testfunction and consider a fixed { G I, when we need to 
inspect the e — > limit of 



(a,W £ [ijt /£ ])= / dp / dA ; o(p 1 A ;| J)e- i ^ w ^*^*- e »)^g(*-ef)^(fc + e|). 

JR3 Jj3 

(4.1) 

First we identify the function 0o which contains the contributions of the long-wave 
excitations. Let 0q be defined by (I3.8p . Since limsup e ^ ||0olU 2 = limsup e ^ HV'olta 
is bounded by Assumption 13. II there exists a subsequence and a function O G L 2 (IR 3 ) 
such that 4>q converges weakly to O m £ 2 (IR 3 ). 

Using a localization function \ which will be specified later we split the rhs of 
(14.1ft into three parts and an error term so that the contributions short-, medium- 
and long-wave excitation can be analyzed separately: 



dp / dA;a(p,fc,f)e- i K-( fc + £ l)--(^f))^(fc- £ |)^(A; + £|)(l- x ), 
I*= I dp [ dk a( P} k, *) e -i!(-( fc +4)— ( fc -4)) (^§(fc _ £ R) - e -§0 o ( * - I)) 



T 3 



x(V>g(* + e|)-e-*to(f + i))x, 
I? = e* ( dp [ dk o(p, fc, 4) e - l K"(*+«f )-"(»-f )) o (f - I) o (f + |) x , 



T 3 

i2 = /-/>- - 7™ 



The definition of implies that 

(a, W^%]> = /> + /? + ir + ^ (4-2) 
To localize the oscillations in Fourier-space we need smooth cutoff-functions. 
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Definition 4.1 Let f e C^(R, [0, 1]) denote a fixed function which is symmetric, 
f(—x) = f(x), strictly monotonically decreasing on [1,2], and 





if \x\ 


1° 


if \x\ 



/(*) 



We define further ip e C (oo) (M 3 , [0, 1]) by <p{k) = f(\k\). 



(4.3) 



Let < p < | be arbitrary and set x{fy — V 9 \Jp) V \p) wnere ^± = k±e^. We 
will continue to use this shorthand notation, under the tacit assumption that k± is 
always really a function of k and ep. All four terms depend on e and p, although we 
do not denote this explicitly in general. 

The first term containing 1 — x in ( 14.21) is zero, if \k±\ < p, while the remainder is 
zero, if \k + \ or \k-\ > 2p. Thus the chosen decomposition splits the integration over 
k and p into "large", "intermediate" and "small" wave numbers. We will demonstrate 
that the following convergences hold: there is a sequence of p and a subsequence of 
e such that 



lim lim i> = 

p^Oe^O 

lim lim = 

p-»0e^0 ^ 

lim lim sup 
lim lim sup \R\ =0. 



dfit(x,k) b(x,k,jq), 



■•xS 2 



dfif(x,q) b(x,0,q), 



Clearly, equations lEOI-IOl imply (1331) . 



(4.4) 

(4.5) 

(4.6) 
(4.7) 



Large wave numbers. 

We split i> further into two parts using 



Mi 



(4.8) 
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Letting h p (k) = 1 — ip y^j , which is a smooth function, the integral then becomes 
J> = Ii, + Ri, where 

Jj= I dp f dk a(p, k, ±)h p (k) e -^("(*+^)-»(k-eZ)) ^ (fc _ £ R) fr Q ( k + e £) . (4.9) 
Jrs Jt 3 

The remainder R\ can be estimated using \(p\ < 1 and 

v(^)-v(j) =± 4 f l dsp.V ¥ >(j(k±s^p)), (4.10) 
which yield the bound, with a universal constant C, 

\R,\ <C|sup||a(-,A;,g)||^ +2 ||^|| 2 ||V^|| 00 . (4.11) 

k,g 

Therefore, there is a constant d such that \R\\ < c'e/p and thus R\ — > when e — > 
for all p. 

We then consider J>. The presence of /i p guarantees that the integrand is zero 
unless \k\ > p. Thus we can change a(p, k,k/e) to b(p, k,k/\k\) in the integrand with 
an error R2 bounded by 

\R2\<C sup \\a(.,k,q)-b(-,k,±)\\ S4+1 \\4\\ 2 (4.12) 

k,\q\>p/e 

where C is a universal constant. Therefore, the assumptions imply that R2 — > 
when s — > for all p. On the other hand, for \k\ > p and \p\ < p/e, inequality flA.2j) 
implies that 

I \ (u(k + el) - cu(k - el)) - p ■ uo{k) \ < C 3 e^ < C^ p \p\ 2 . (4.13) 

Therefore, using the estimate |e lx — e iy | < min(|x — y|,2), valid for all x,y G M, 
we find that we can further change the t-dependent exponential in the integrand to 
e -itp-Vu(k) w j^j 1 a error r 3 bounded by 

\R 3 \<C sup \\b(;k,q)\\ S4+3 (f\t\+ [ dp(p)- d - 3 ) (4.14) 

k,q V J\p\>p/e / 

for some constant C. Therefore, also lim e ^ R3 — for all p. In summary, J> = 
+ _R 2 + -R3, where i? 2 , -R3 are negligible, and 

i*=i*(e,p) = [ dp [ dk b(p, k, *\h p {k)e- lt v^ r (k - ej) ft (k + e^ 2 ) . 
Jr 3 Jt 3 

(4.15) 
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Let us for a moment consider the lattice Wigner transform Wf att of i/jq, as defined 
in [Tj. As pointed out after Definition 12.51 then for any testfunction / G S(M? x T 3 ), 
we get an admissible test-function by the formula a/(x, k, q) = f(x, k) and then also 
(/) ^latt) = ( a f^ £ bPo\) ■ Since ipQ is a norm-bounded sequence, the sequence W( &tt 
is weak-* bounded, and thus there is W^ tt G iS'(M 3 x T 3 ) and a subsequence along 



which WU A < tt . 



Since the sequence ip^ is by assumption also tight on the scale e -1 , we can then 
apply Theorems B.4 and B.5 of [7] and conclude that is given by a positive, 

bounded Radon measure /i on M 3 x T 3 such that for all continuous functions / G 
C(T 3 ) and p G M 3 , 

lim/ dfc/(fc)^g(fc-e§)^(fc + e§) = / dp(x, fc) f (k)e~ 27ripx . (4.16) 



As Va;(A;) and fc/|A;| are continuous apart from k = 0, the function /c i— > h p (k)b(p, k, jjMe l4 P -Va; ( fe ) 
is everywhere continuous for all p > and p G M 3 . Therefore, by the dominated con- 
vergence theorem, for all p we find 

limij= / dp / d/i(x, fc) S(p, k, ^Ah p (k)e-^ ip < x+tv ^ 2 ^ 

= / d/i(x, fc) h p (k)b(x + t±Vuj{k), k, g . (4.17) 

jR3xT3 

When p — > 0, the integrand approaches pointwise b(x + t^Va;(A;), fc, m) apart from 
= 0, when the limit is 0. Therefore, by the dominated convergence theorem 

lim / d{i(x,k)h p (k)b(x + t±Vuj(k),k,^) 

= [ dn(x,k)b(x + t±Vu(k),k±) = [ d// t (x,fc)6(x,fc,|i). (4.18) 

JR3 X T3 ' ' JlR3 xT 3 1 1 

where we have defined the bounded, positive Radon measure fi t using p = p|r3 x t3 
in the formula (13.51) . We have shown that equation f)4.4p holds. 

Small wave-numbers and the remainder 

After a change of variables q = - one obtains that 



'eq- 



I- v = / dp / dga(p,^,g)e- i ^^+)-^-))0 o (g-)0o(g + )^( £ ^)^(- 

(4.19) 
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where q± = q ± |. We can immediately replace the integration region for the q- 
integral by M 3 . To see this, note that the integrand is zero, unless \q ± || < ^ for 
both signs. Since 2q = q + + g_, this can happen only if also \q\ < — < ^, which 
implies that the integrand is zero if ||g||oo > 
However, if \eq\ < 2p, then 

eq, q) - a(p, 0, g) | < sup | V*o(p, fc, g) | 2p, (4.20) 

and thus we can replace in the integrand the function a(p, eq, q) by a(p, 0, q), with an 
error R[ which is bounded by Cp with a constant C independent of e and p. Thus 
we only need to consider the integral 

Ji= / dp / dg a(p, 0, g) e-'K^*)""^)) F^(g_)F £ ^(g + ), (4.21) 

where F £,p {q) = 0o(g)<p(^)- Next we use estimate (1A.3j) . which implies that, if now 
\p\ < i, then 

|i(w(eg + )-w(eg_))-wo(g + )+wo(?-)| < C 4 e|p| |g| . (4.22) 

Following the same argument as earlier, we can then conclude that the t-dependent 
exponential can be changed to e _1 *^°^ + ^ wo ^ _ ^, with an error R' 2 which satisfies the 
estimate 

|i^|<Csup||a(-,0,g)|| S)d+2 (|t|p+ / dp(p)- d - 2 ). (4.23) 

Thus limp^o limsup e ^ \R' 2 \ = for all t. Finally, we need to change F £,p to 4> , with 
an error R' 3 which can be bounded by C\\F e,p — <p \\. Since the bound goes to zero 
when e — > and 

IZ V = [ dp [ dq a(p, 0, g) e ~^ + )-^-)) o ( g _)0 o(g+ ) + R[ + R ' 2 + ( 4 .24) 

equation ( 14. 6ft has been established. 

Similar estimates can be employed to demonstrate the vanishing of the remainder, 
equation (14 .70 . From the definition of R we get 



R= dp dqa(p,£q,q)e- [ ^^-^~» feg_) - o (g_)) o (g+) 

JR 3 </T 3 /e V 7 

X <p(^±) tp(?±) + f dp f dqa(p,eq,q)e- [ ^ u){£q+) - uj{eq - )) 

x Uql(k(q + ) - 0o(g + )) <fi( £ -f) (4.25) 
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We then apply the above estimates to remove the e-dependence from all other terms 
in the integrands, apart from the differences 0q — 0o- The error has a bound which 
vanishes when p — » 0. We are then left with 

[ d v [ dCaiv-ZAUv + O)^ ^" ^ (k(0 - MO) Mv) 

JM 3 J R 3 V / 

+ f drj [ d£ a(ri - I 0, |(r/ + 0) e-'^W-" ®) J^)U £ (v) ~ Mvj) , (4-26) 
Jm. 3 Jm. 3 v / 

which vanishes as e — > 0, since 0q converges weakly to <f>Q. This establishes (14.71) . 

Intermediate wave-numbers 

Changing coordinates q — ~ yields that 

J== / dp / d? o(p, eq, q) e -f (<"N+)-<^-)) f^(q_)f^( q+ ) , 
Jm 3 Jt 3 /e 

where f £ ' p (q) = (p (0g(g) - 0o(g)). Let M > be arbitrary. We split off the 

values \p\ > M from the integral defining /<. The difference R' A = R' 4 (e, p, M) can 
be bounded by C f^ >M dp (p)~ d ~ 1 and thus liniAf^oo sup p £ \R' A \ = 0. We divide the 
remaining integral over \p\ < M further into two parts using the identity 

l = 0(£)*(&)+l-0(fe)0(&). ( 4 - 27 ) 

where (p — 1 — (p. If |g| > 5M, then \q±\ > 4M and the second part is zero. It can 
be checked by inspection that the sequence (f £ ' p ) £ is bounded and tight, and it has a 
weak limit zero. Of these properties only the tightness is non-obvious, but this can 
also be easily deduced from the formula 

f E ' p (x) = pV 3 / 2 V r o (i)0(p(i - H) - / d yM* + ~y)0(y)- (4.28) 

1& 3 " J* 3 

Therefore, by Lemma IA.21 lim e ^ \\f £ ' p \\L 2 (B 6M ) = f° r a ^ M,p. This implies that 
the contribution of the second term, denoted by R' 5 , satisfies lim e ^ |i? 5 | = for all 
M,p. 

We are thus a left with 

I*= f dp [ dqa{p, 0, q) e --<Mi+)-M<!-)) g^ M (q_)g E 'P' M (q + ), (4.29) 

J\p\<M JM 3 
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where g £,p,M (q) = <p (^) f £ ' p (q) and the integrand can be non-zero only for |g| > M. 
We thus only need to consider \q\ > M and \p\ < M. First we replace in the integrand 
0, q) by b(p, 0,q), q = A, with an error i? 6 which is bounded by 

|i? 6 |<C sup \\a(;k,q)-b(;k,M\ S4+1 (4.30) 

fc,|?|>M 

for some constant C. Then we change er lt( ^^i+^~ w °ii-^) to e~* ltp " Va ' W with an error, 
R' 7 which can be estimated using ( 1A.4[) which proves that there is a constant C such 
that 

1^1 <C&suv\\b(;k,q)\\ s>d+3 . (4.31) 

fc,q 

Therefore, J* = I% + R' 6 + R' 7 , where 

7j= / dp I dqb{p, 0, g) e - itp - Ww °^ gt<P< M {q„)g £ ' f) ' M {q + ). (4.32) 

and liniM-nx) su P p , e |-R 6 + R' 7 \ = for all t. 
Let 

fetfo g) = b(x + t£Vw ($), 0, g). (4.33) 

Then b t G iS(M 3 x S 2 ), and it has an extension to a function J t G <5>6, i- e -> there is Jt 
such that Jf(x, q) = b t (x, q) for all |g| = 1. On the other hand, then 

I%= f dp f dqb t (p,q)g^(q_)g £ > p > M (q + ) = (J t ,A £ ' p ' M ) (4.34) 

J\p\<M JU 3 

where A e ' p ' M G S' 6 denotes the distribution 

J„(J i A<>»")= [ dp I dqJ(p,q)g £ ^(q^)g £ ^ M (q + ). (4.35) 
J\p\<m Jm. 3 

Clearly, each A £,P,M has support in R 3 x S 2 , and there is a constant C such that 
for all J,e,p,M 

|<J,A £ '^}| <C\\g £ ' p ' M \\ 2 \\J\\s, d+ i- (4-36) 

However, since we have ||g e ' p ' M || < ||/ £ ' p ||, where ||/ e ' p || is bounded in e, Banach- 
Alaoglu theorem implies that the family (A E ' p,M ) EtPt M belongs to a weak-* sequentially 
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compact set. Therefore, for every p, M there is a subsequence of (e) and A P,M such 
that A £ ' P ' M — A P,M along this subsequence. In addition, for every p there is A p and 
a sequence of integers M such that A P,M — A p along this sequence. Finally, there is 
A e 5g and a sequence of integers iV > 4 such that for p = jj, A p — ^ A. 

All of the above distributions clearly must have support on IR 3 x S 2 . We will soon 
prove that, in addition, for all / G S 6 and p 



{\f\ 2 ,A p )>0. 



(4.37) 



This implies then that also (|/| 2 ,A) > 0. Therefore, by the Bochner-Schwartz the- 
orem, there is a positive Radon measure /i H on IR 6 such that for all testfunctions 
J, (J, A) = f p u (dx,dk) J(x,k). Since also /i H must have support on IR 3 x S 2 , we 
can thus identify it with a positive Radon measure //q on I 3 x 5 2 . By considering 
testfunctions J(x, k) = e~ s x in the limit 5 — > 0, it is also clear that /i^ must be 
bounded. We then define the positive, bounded Radon measures pf , t G R, by the 
formula (13.61) . It follows from the construction of /x^ that equation (14.51) holds along 
the above sequences p, e. 

The main missing ingredient is provided by the following Lemma 

Lemma 4.2 Forp,q G M 3 , Zet denote q± = q ± |, g± = g ± /|g_|-| ; and g = g/|g|. 
There is a constant C such that for all q, g' G I 3 , and f E Sq 



dxe 



-2-Kip-x 



■ x f(x,q')f(x,q) 



<C(p) 



-d-l 



2 

S,d+V 



(4.38) 



If, in addition, \q\ > M and \p\ < M, then also 



^ p (\f\ 2 )(p,q)- / dxe- 2 ^ x f(x,q + )f(x,q.) 



" "IW (4-39) 



Before proving the lemma we demonstrate that it implies inequality ([4.370 . Let then 
/ G Sq be arbitrary. Define g±, g± as in Lemma I4~2l except that let here also = 0, 
and let 



je,p,M,f 



dp / dq 



dx e 



f(x,q + )f(x,q_) g^ M (q_)g^ M (q + ). (4.40) 



Note that this integral is well-defined by (14.381) . Then, by estimate (14.391) and uniform 
boundedness of g e ' p ' M , there is a constant c such that 



|(|/| 2 ,A £ ' P ' M ) -p,pWJ\ < 



S,d+2 



dp(p) 



\p\>M 



-d-l , l_ 



(4.4i; 
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On the other hand, always I £ 'P' M 'J > 0. To see this, consider first the case when 
ge,p,R £ g Then g £ 'P' R g S, and by changing variables from (q,p) to in 
(14.401) . and then using Fubini's theorem to reorder the integrals, we find that 

F' p ' M ' f = [ dx\G(x)\ 2 , withG(x)= I dqe 2niq x f(x,q)g £ ' p ' R (q) E L 2 . (4.42) 

Since I £ >p> M 'f depends L 2 -continuously on g £ 'p' R this implies that also for general 
g£,p,R ^ we have I E >P> M >f > 0. Since the right hand side of (14.411) vanishes if first 
e — > and then M — > oo, we must thus also have (|/| 2 , A p ) > 0. This proves ( 14.371) . 

The only remaining task is to prove Lemma [4.21 Consider first (14.381) . If \p\ < 1, 
we have trivially a bound 

/ dx\f(x,q')\ \f(x,q)\ < \\f\\l d+1 I dx(x)- 2d - 2 <c\\f\\l d+1 . (4.43) 

If \p\ > 1, we perform N partial integrations in the direction of p, that is in the 
direction p = p , yielding 

/ dxe- 2 ^JM)f(x, q) = [ dxe~ 2 ^{p ■ V) N (lM)f(x, q)) . 

(4.44) 

By the Leibniz rule, 

(p " V)" (7M)f(x, ?)) = X) ( ) (P • V) n TM)(P ■ V) N ' n f(x, q) (4.45) 

n=0 ^ ' 

which is bounded by c' N (x)~ N \\f\\s^ N . Choosing N = d + 1 then yields (14.381) for 
some constant. Adjusting the constant C so that the bound is true also for \p\ < 1 
proves that (14.381) is valid. 

To prove ( 14.391) . consider q, p as required in the Lemma. Let q(s) = q + |p, when 
\q(s)\ > \q\/2 > for all \s\ < 1, and q(s) is thus well-defined and smooth. Therefore, 
for any g G Sq and So G [—1, 1], 

/•so 

g(x, q(s )) - g{x, q) = / ds £g(x, q(s)) 

Jo 

= / ds {W\P - ■ V ^> q) \ q=4(s) (4.46) 

J 
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implying 

\g(x,q ± )-g(x,q)\ < 2§ sup | V 9 g(x, q)\. (4.47) 

l?l=i 

Also for all g±, g 2 e Sq , 

gi(x, q)g 2 (x, q) - gi(x, q-)g 2 (x, q+) 

= (gi(x, q) - gi(x, q-))g 2 (x, q) + #i(x, g_) (g 2 (x, q) - g 2 (x, q + )) . (4.48) 

Following the steps made in the first part of the proof, and replacing the earlier 
estimates with the above more accurate ones when necessary, we can conclude that 
the constant C can be adjusted so that for these values of q,p also (14.391) holds. This 
completes the proof of equation (14.51) . 



Energy equality 

The energy equality (13. 71) follows by considering a sequence of testfunctions a s = 
e -8 x anc j t a ]d n g 5 — > 0. To see this, first note that the right hand side of ( 13. 71) 
is clearly independent of t, and thus it is enough to consider t = 0. Thanks to 
equation (I2.26P and to the tightness of the sequence ipQ we obtain for this particular 
testfunction that 

hm hW, W e [^]> = hmhm(aV) = lim||^|| 2 

<5— >0 £^0 5^0 e— +0 e^O 

Equation (13.31) implies that 

\im(a s ,W £ [ij £ }) = [ dxa s \Mx)\ 2 + I [ a s (x) dfi (x, k) + [ [ a\x) d${x, k). 

JR3 J R 3 J T 3 J R 3 J S 2 

Sending 5 to yields that 

lim \\r \\ 2 = \\Mh m + M^ 3 x T 3 ) + /,«(R 3 x S 2 ), 
and the energy equality has been established. This finishes the proof of Theorem 13.21 



5 Proof of Corollary 13.4 

Let Iq denote the original sequence of e, and consider an arbitrary subsequence I of Iq. 
Since (ipQ) ee i then also satisfies Assumption 13. H we can conclude from Theorem 13.21 
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that for every / there is a subsequence I' such that ( 13.31) holds for all t with the 
initial conditions given by some triplet (fii, [if, <fri). From the construction of the 
subsequence in the proof of Theorem 13.21 we know that <pj can be chosen as the weak 
limit of 4>q along the subsequence I'. The first assumption thus implies that we can 
always choose (pj = O - Let us also denote fi = fij and /1q = fif Q , and to prove the 
stated uniqueness, we will prove that /// = /j,q and [if = $ for all /. 

Consider first an arbitrary a G Cc^^M 3 xT 3 ). Let a(x, k, q) = a(x, k) for k 7^ and 
define a(x, 0, q) = 0. Then a is an admissible test-function with ao = = b(k, 0, q) 
and for any subsequence I we thus obtain, using the second assumption 



/ dfii(x, k) a(x, k) = \im{a,W e [ipQ}) = / dfi (x,k)a(x,k). (5.1) 

jR3 xT 3 SGI ' - jR3 xT 3 

Such a are dense in C C (1R 3 x T 3 ), and thus //j = /j,q on IR 3 x T 3 . 

Consider then an arbitrary b G Ci (K. x S 2 ). Let <y9 be a smooth cutoff function 
as in Definition 14.11 and define a(x,k,q) = b(x,q/\q\) (1 — tp(2q)). Then a is an 
admissible testfunction and indeed liniR^oo a(x, k, Rq) = b(x,q) for all |g| = 1. By 
the already proven results we then find that for any subsequence I 

/ dfif(x,q)b(x,q) = hm(a,W e [ip E ]) - / dn (x, k) b(x, A) 

JR3 XS 2 e^O JR3 X T3 

- (a ,w£iM) = [ difi(x,q)bfcq). (5.2) 

JRZxS 2 

Therefore, also fif = which concludes the uniqueness part of the proof of the 
Corollary. 

Finally, define for t 7^ the triplet (// t , fif, 4> t ) using (fi , //q , <f> ) as the initial 
data. By the above proved uniqueness, for any subsequence /, (13.31) holds along the 
subsubsequence I'. As the right hand side is thus independent of J, this proves that 
the limit also holds along the original sequence Jo- This completes the proof of the 
Corollary. 



A Technicalities 

The proof of Theorem 13.21 relies on two simple lemmas which are provided here. The 
first lemma summarizes several properties of regular acoustic dispersion relations. 

Lemma A.l Let u be a regular acoustic dispersion relation and recall the definitions 
of X, Aq and Uq. Then the following assertions are true. 
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1. cue C( 3 )(T 3 \{0},M). 

2. VA(O) = and A > 0. 

3. There are constants C\,C2 > such that for all \k\oo < |, 

w(ife) > Ci|&|, |VA(fc)| < C 2 |fc| . (A.l) 
In addition, \\Vu)\\oo < oo. 
4- There is C 3 such that, if e > 0, p G R 3 ; and | A;|oo < |, wnt/i |fc| > e\p\, then 

\cu(k + iep) - uj(k - |ep) - sp ■ Vuj(k)\ < C^e 2 ^. (A.2) 

5. There is C4 such that, if e > and 6 R 3 with \p\, \\q\\oo < 2 £_1 > ^ en / or 

|w(eg + ) - u(eqJ) - w (eg + ) + LO (eq_)\ < C 4 e 2 \p\ \q\ . (A.3) 

6. There is C§ such that, ifp, q G R 3 wift g 7^ and \p\ < \q\, then for q± = q±\p, 

-wb(g_) - p- Va; (|J)| < C 5 ^. (A.4) 

Proof: The first item is obvious, and the second one follows from the assumptions, 
since is a minimum. The second inequality in flA.lj) follows by using item [2] and 



||-D 2 A||oo < 00, when by Taylor expansion |VA(fc)| < C 2 \k\ for all k G R. To prove 
the first inequality, we first note that, by continuity, also ||D 3 A|| 00 < 00. Thus there 
is d such that for all k G R 3 



|A(fc)-A (fc)| <c'|A;| 3 , (A.5) 

where Aq(&) = \k ■ A k. Since A Q > 0, there is c > such that Xo(k) > c 2 k 2 for all 
k. Thus there is 5' > such that for all \k\ < 5', we have |1 — X(k)/Xo(k)\ < |, and 
for these k therefore also u(k) = \/X{k) — Xo(k) + Xo(k) > |yAo(A?) > §|fc|. Since 
A has no zeroes in the complement set, the constant can then be adjusted so that 
(|A~TD holds for all k with \k\oo < §. 

We still need to prove the third property, boundedness of Vu. For later use, let 
us, more generally, consider a non-neg ative / G C (2) (R 3 ), q G R 3 , and k such that 
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f(k) 0. Then we have 



(q ■ V) yW) = ^7y^y<? * V ^)> ( A - 6 ) 
(q • V) 2 ^ = "4^(9 • V/W) 2 + ^==(g • V) 2 /(A0- (A.7) 
This implies that there is a constant C such that for all g e K 3 and k 7^ 0, with 



\k\ < 2- 



\q ■ Vuj(k)\, \q ■ Vu (k)\ < C\q\ , (A.8) 
|(g-V) 2 ^(A;)|,|(g-V) 2 a;o(A:)|<C^. (A.9) 



In particular, by periodicity therefore ||Vu;| 

To prove item HI consider k,p,e as in the claim, and define a function f(s) = 
u)(k + ) — u(k_) with k± = k ± s^p and s G [—£,£]. Then \k±\oo < §|fc|oo and 
\k±\ > \k\ ~ e \v > \\k\ > 0, and thus / belongs to C^ 3 >. In particular, /(0) = 0, 
f (0) =p. Vw(fc), and 

™ 4 [^b (I' " ^) (f ■ V ) V) 

? • VA(A; + )) 2 + ^-i-^ • VA(M) 2 j • (A.10) 



2w(A; + ) 3 V2 v +, J 2w(A;_ 

Using item[3], we find that there is C such that this is uniformly bounded by C|p| 2 /|/c|. 
Then a Taylor expansion at the origin proves item HI 
Since for all k ^ 0, 



V A(fc) + V Ao(fc) 



(A.11; 



and then also 



q- V(u(k)-u (k)) 

_ q ■ VA(fc) - q ■ VAp(fc) g- Vu;(A;) + g ■ V^ (A;) 

u;(fc) + u; (k;) (u;(fc) + u; (fc)) 

By Taylor expansion at the origin, we find that there is C such that for all q G M 3 , 

|A(fc) - A (fc)| < C"|A;| 3 , \q-V\{k)-q-V\ {k)\<C'\k\ 2 \q\, (A.13) 
|(g • V) 2 A(fc) - (q • V) 2 A (fc)| < C\k\ \q\ 2 . (A.14) 
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Thus there is also C such that for all q £ M 3 , and \k\oo < f , with k ^ 0, 

\q-V{co(k)-cu (k))\<C\k\\q\, \(q ■ V) 2 (u(k) - u (k))\ < C\q\ 2 . (A.15) 

Let us then consider q, p, e satisfying the assumptions made in the final item. Since 
then ||£<7±||oo < §, we can apply the previous estimates. If p = 0, then g_ = q + and 
the bound in flA.3j) is trivially valid for any C4. Consider thus p 7^ 0, and assume first 
that p is not proportional to q. Since then the line segment [0, 1] 3 s 1— > eg + sep/2 
does not pass through the origin, the function s 1— > ui{eq + sep/2) — uo(eq + sep/2) 
is in C^([0, 1]). We make a Taylor expansion of this function at s — 0, yielding 

cu(eg + e|) - uj (eq + e|) = u;(eg) - ^(eg) + e f ' (^ 7ll, ( £ '?) _ Vcj (eg)) + -R . 

(A.16) 

Here ( 1A.15I) implies \R\ < Ce 2 \p\ 2 , since 

ds (1 - s) (el ■ v) (cu(eg + sef ) - cj (£? + sef)) . (A.17) 

This proves that (1A.3j) holds in this case for some constant C4. In final remaining 
case p oc q, we choose a direction it orthogonal to g, and use the previous estimate 
with p + 5u instead of p for an arbitrary < 5 < 1. Since the left hand side of f[Q|) 
is continuous in 5, the bound must then hold also for 5 = 0, proving the validity of 
the estimate also in this case. 

To prove (1A.4j) . we use the fact that by assumption \q±\ > ||g| > 0, and thus 
denoting q = q/\q\, we get 

( V P \ 

u (q+) - wo(g-) = |g| I w (g + ^j) - ^ (g - ^j) J =p • V^ (g) + R (A.18) 

where < C|p| 2 /|g|. We have thus completed the proof of the Lemma. □ 

The second lemma recalls a well-known fact in Fourier analysis: the Fourier transform 
of a bounded and tight sequence of L 2 functions converges strongly on compact sets. 
For the convenience of the reader we give a proof. 

Lemma A. 2 Let f e e L 2 (IR d ) be a bounded and tight sequence of functions such 
that f £ ^0 as e — > 0. Then 

Um 11/1^(0) =0. (A.19) 
/or every Qcl d with a finite measure. 
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Proof: Choose an arbitrary number (3 > 0. By tightness of f £ there exists a number 
Rp > such that 



limsup / dx\f £ (x)\ 2 < (3 2 . (A.20) 

e^O J\x\>R!3 



Define now 



fl^ f' 1 ' " W " % (A.21) 
I else. 

By the boundedness of f 6 ' 13 in L 2 (IR d ) there exists / and a subsequence (V) such 
that f 6 '' 13 g 13 in L 2 (R d ) as e' -> 0. Estimate ([QOjl implies that limsup e ^ ||/ e - 
/ e '^lli 2 (R 3 ) — /^ 2 - Since also f £ — 0, we have 

||/|| 2 = lim|(/,f'^>| < ll/H limsup Hr-rt (A.22) 

and, therefore, || 1 1 z, 2 (M d ) < /5- 

Since g 13 has support in a ball of radius Rp, also gf (x) = —2mxg^{x) G L 2 (M. d ). 
Therefore, g 13 G /^(R 4 ), and Vg 13 = g{. Similarly, also G if 1 ^) and it is 
straightforward to check that V/ £ — V/ in L 2 . Since f2 is assumed to be a 
set of finite measure, we have then lim e /^ ||/ £,/3 — ,r IU 2 (n) — (for a proof, see 
for instance Theorem 8.6. in [15]). This result yields the following estimate for the 
original sequence f £ : 

limsup \\f £ '\\ L 2(n) 

< limsup ||f ' - + limsup ||f ' ? - f?\\ L , {n) + \\f%2 (m < 2(3. 

£'^0 £'^0 

Since (3 can be arbitrarily small, we obtain that there must be a subsequence e" such 
that linv/_> ||/ e "|U 2 (n) = 0. 

Since the assumptions on the sequence f £ are preserved for subsequences, we can 
consider an arbitrary subsequence and apply the above result to it. Then we can con- 
clude that for every subsequence there is a subsubsequence along which ( 1A.19I) holds. 
This implies that the limit ( 1A.19I) actually holds also along the original sequence. □ 
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